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COThe NF -kB signalling system is involved in a 
•^firiety of cellular processes including immune re- 
'wponse, inflammation, and apoptosis. Recent ex- 
periments have found oscillations in the nuclear- 
i-eytoplasmic translocation of the NF-kB transcrip- 
mjni factor [Hoffmann et al. (2002) Science 298, 1241; 
^Ison et al. (2004) Science 306, 704.] How the cell 
j*^es the oscillations to differentiate input conditions 
•artd send specific signals to downstream genes is an 
'TJjpen problem. We shed light on this issue by exam- 
ining the small core network driving the oscillations, 
L wnich, we show, is designed to produce periodic 
T agikes in nuclear NF-kB concentration. The oscil- 
lations can be used to regulate downstream genes in 
ta""variety of ways. In particular, we show that genes 
id whose operator sites NF-kB binds and dissociates 
jKst can respond very sensitively to changes in the 
fmput signal, with effective Hill coefficients in excess 
Uff)20. 

o ■ 

— NF-kB is a family of dimeric transcription factors which 
.Qrticipatc in the regulation of a number of cellular pro- 
r£ekses including immune response, inflammation and apop- 
^}s m^UHlEl- Extensive experiments using electrophoretic 
mobility shift assay and single-cell flourescence imaging have 
.^Jind oscillations in the nuclear-cytoplasmic translocation 
^Jthe NF-kB transcription factor in mammalian cells |31|7Jj, 
fcith a time period of the order of hours. NF-kB can be 
, activated by a number of external stimuli [Jj including bac- 
teria, viruses and various stresses and proteins (e.g., tu- 
mor necrosis factor-a, TNF-a, which was the signal used in 
refs. In response to these signals it targets over 150 

genes including many chemokines, immunorcccptors, stress 
rcponse genes, as well as acute phase inflammation response 
proteins [Zj. Experiments show that NF-kB docs not regu- 
late all its downstream genes in the same way. For example, 
the chemokine gene RANTES turns on much later than an- 
other chemokine IP-10 after TNF-a activation Thus, 
the two main questions raised by the dynamics of the NF- 
kB system are: how does the network of interactions pro- 
duce oscillations, and how does the cell use the oscillations 
to differentiate input conditions and send specific signals to 
downstream genes? In this paper, we elucidate the small 
core network driving the oscillations and show that it is de- 
signed to produce periodic spikes in nuclear NF-kB concen- 
tration. We show that the spiky oscillations are extremely 



robust to variation of parameters. We further argue that 
the spikiness is associated with an increased sensitivity of 
the system that could be used for differentially regulating 
downstream genes. 



Extracting the Core Feedback Loop 



Hoffman et al. have constructed a long list of chemical reac- 
tions between 26 different molecules in the NF-kB system, 
including reaction constants We reduced this system 
to the core feedback loop (Fig. Qj3) generating oscillations. 
The reduction was done in three steps: the first, removing 
molecules which have no feedback from NF-kB and deleting 
slow reactions where faster alternate pathways exist (e.g., 
export of nuclear NF-kB) resulted in a 7-variable model. 

The interactions in this model arc schematically displayed 
in Fig. It consists of cytoplasmic and nuclear NF-kB, 
its inhibitor, IkB, and IkB kinase (IKK) which phosphory- 
lates the inhibitor, leading to its degradation. The inhibitor 
forms a complex with NF-kB which, in the cytoplasm, pre- 
vents its transport into the nucleus. Only free nuclear NF- 
kB is imported into the nucleus. In contrast, from inside 
the nucleus, only the complex can be exported, not the free 
NF-kB. IkB is known to occur in several isoforms. Cells con- 
taining only the IftBa isoform show sustained oscillations, 
while cells with only the IkB/3 or e isoforms do not show 
oscillations. Wild type cells, with all three isoforms, typi- 
cally exhibit damped oscillations [S] . The difference between 
these isoforms is that only InBa is activated by NF-kB 015]. 
In contrast IkB/3, e are produced at a rate independent of 
NF-kB and so lie outside the feedback loop (see support- 
ing text for the equations governing the dynamics of the 
7-variable model.) 

Coarse graining over fast chemical reactions involving 
complex formation reduced this system to 4 variables. Fi- 
nally, based on numerical observations, we found we could 
effectively eliminate nuclear IkB, giving the model in Fig. 
03. More details of the reduction process are given in sup- 
porting text. 
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Figure 1: Schematic diagram of interactions in the NF- 
kB signalling system. A. 7-variable model. B. 3-variable 
model. 

Results 

3- Variable Model of NF-kB Oscillations 

The core feedback loop, we find, consists of only three con- 
stituents (Fig. IB): nuclear NF-kB (N n ), cytoplasmic IkB 
(I) and IkB mRNA (I m ). NF-kB dimers activate produc- 
tion of IkB mRNA which translated to IkB inhibits nuclear 
NF-kB production, completing the feedback loop. 

The dynamics of the system in Fig. ^3 is captured by 
three coupled ordinary differential equations: 
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For ease of analysis we have rescaled all variables to be di- 
mensionless (the original equations and the rescaling pro- 
cess are described in supporting text.) A, B, C, 6, e are di- 
mcnsionless parameters dependent on the reaction constants 
(see supporting text). The external signal is supplied by 
IKK that enters the equations through the parameter, C, 
which is proportional to IKK concentration. The model has 
no spatial degrees of freedom. It has been suggested that 
flow induced stresses can trigger NF-kB mediated gene ac- 



Figure 2: Sustained oscillations of nuclear NF-kB (N n ), red, 
and cytoplasmic IkB, green, for A = 0.007, B = 954.5, C = 
0.035, S = 0.029 and e ~ 2 x 10~ 5 (using parameter values 
taken from jS], see supporting text). The shape, phase and 
period are remarkably similar to the experimental plot in 
Fig. 3B, left panel, in ref. [H]. 



tivity JO] but here we assume that each compartment, the 
cytoplasm and nucleus, is well-mixed. 

The first term in equation (la) represents the import of 
free cytoplasmic NF-kB (whose concentration is 1— N n ) into 
the nucleus. This is hindered by the presence of cytoplasmic 
IkB which sequesters NF-kB in the cytoplasm. Parameter 
A is proportional to the NF-kB nuclear import rate. The 
second term in (la) derives from the export of NF-kB from 
the nucleus via the NF-kB-IkB complex, which is why the 
term also depends on /. B is proportional to IkB nuclear 
import rate. <5 sets the concentration at which half the nu- 
clear IkB is complexed to NF-kB and it depends both on 
the rates of association and dissociation of the complex as 
well as the export rate. 

The first term in equation (lb), the rate of production 
of IkB mRNA, contains the square of N n because the pro- 
duction is activated by NF-kB dimers^. The second term is 
the degradation of the mRNA whose rate, in these rescaled 
equations, sets the overall timcscalc. It is easy to modify 
this equation to deal with the j3 and e isoforms of IkB sim- 
ply by adding a constant for their NF-KB-independent rate 
of production. 

Equation (lc) has one term for the production of cytoplas- 
mic IkB from its mRNA and a second for its degradation 
due to the presence of IKK. This degradation is proportional 
to the concentration of the NF-kB-IkB complex, which de- 
pends on both / and (1 — N n ), the concentration of cyto- 
plasmic NF-kB. e sets the concentration at which half of the 
cytoplasmic NF-kB is in the complex. C is proportional to 
the rate of degradation and to the IKK concentration. 

Fig. [2] shows a plot of nuclear NF-kB and cytoplas- 
mic IkB concentrations obtained in simulations, using 
parameter values from ref. 0. Our model predicts the 
following experimentally observed facts El : i) sustained 
oscillations in cells with only the a isoform of IkB, ii) 



TRef. 1111 argues that making this term linear in N n is better. 
We have checked that this change does not alter our conclusions (see 
supporting text.) 



damped oscillations in wild type cells which include other 
isoforms of IkB, iii) time period of the order of hours, iv) 
spikiness of nuclear NF-kB and asymmetry of cytoplasmic 
IkB oscillations, v) phase difference between NF-kB and 
IkB, vi) lower frequency upon increased transcription of 
IkB. 



Saturated Degradation of IkB is Crucial for 
Oscillations 

A key element in our 3- variable model is the saturated degra- 
dation of cytoplasmic IkB in the presence of IKK (second 
term in eq. lc) due to the Michaelis-Menten complex forma- 
tion between NF-kB and IkB - a complex needed for IKK 
triggered degradation of IkB. The same complex inhibits nu- 
clear NF-kB production because only free cytoplasmic NF- 
kB is imported into the nucleus. 

A stability analysis of the system shows the importance 
of the saturated degradation for oscillations. In general, the 
system has a single fixed point where all concentrations are 
unchanging in time. For small values of e, corresponding 
to strong saturation of the degradation, this fixed point is 
unstable and the system goes into a periodic cycle. As e is 
increased the fixed point becomes stable and the oscillations 
disappear (Fig. |3J). This happens when the value of e be- 
comes comparable to the steady state value of /, which is 
precisely when the degradation rate stops being saturated. 
The saturation is crucial for oscillations because it puts an 
upper limit to the degradation rate, allowing IkB to accu- 
mulate and stay around longer than with the more usual 
/-proportional degradation rate. This effectively introduces 
a time delay into the feedback loop. Negative feedback with 
time delay is known to easily produce oscillations, and this 
mechanism has been used to model the p53 ^2] an d Hes 
oscillations [T^]. Here, instead of being put in by hand, the 
time delay arises more naturally through the dynamics of 
the system. 

Interestingly, the NF-kB core in Fig. ^3 is similar to an 
early model showing oscillations by negative feedback |14| 
which introduced precisely the same kind of saturated degra- 
dation to mitigate the unreasonably large Hill coefficient in 
an even earlier model of Goodwin |15| . The importance of 
this mechanism has also been recognized by Goldbeter who 
has used it in models of various cellular oscillations, e.g., 
the cell cycle [T2|, development in myxobacteria ^7], yeast 
stress response QH] an d the mammalian circadian clock . 
Saturated degradation has also been implicated in models of 
calcium oscillations in cells [201 We further note the p53 
system, which also shows oscillations, contains the degrada- 
tion of p53 via the formation of a complex with Mdm2 [22] , 
which could result in saturated degradation. This suggests 
that saturated degradation might be a very general mecha- 
nism, easily implemented by complex formation and used by 
cellular processes to introduce time delays where necessary. 





unstable 






"T 


stable 




e = 5x 1(T 4 








e= 1.5x 10 " 3 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 




0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 

N n (dimensionless) N n (dimensionless) 



Figure 3: Importance of saturated degradation of /. A. 
Steady state solution of equation (1) for / as a function of 
the parameter e. Dashed line shows the region where the 
fixed point is unstable and the solid line shows the region 
where it is stable (see supporting text). Black dots mark 
the parameter values used in (B) and (C). The crossover 
from unstable to stable occurs where e becomes comparable 
to the steady state value of /, i.e., where the degradation 
of / stops being saturated. B. Trajectory of the system in 
the I-N n phase plane, for e = 5 x 10~ 4 , which converges 
to a stable limit cycle. C. Trajectory of the system, for 
e = 1.5 x 10 -3 , which converges to a stable fixed point. 

Spiky Oscillations and Control of Down- 
stream Genes 

Robust spiky oscillations. One feature of our model 
is that it can produce sharp spikes in nuclear NF-kB. 
This is unusual for oscillations driven by negative feedback 
|23| . We quantify spikiness using the following measure: 
Z = (max(A„) — min(iV„))/mean(iV n ). Oscillations with 
Z > 2 we term spiky, and oscillations with Z < 2, soft. Fig. 
0] shows an example of each type of oscillation, generated 
from the 3-variablc model using different parameter values. 
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Figure 4: Spiky and soft oscillations. Plots of nuclear NF- 
kB concentration as a function of time from equation (1) 
with A,B,C and S the same as in Fig. [21 A. Spiky oscil- 
lations with e = 2 x 10 -5 , the same as in Fig. [21 B. Soft 
oscillations with e ~ 5 x 10 -4 . 
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Figure 5: Robustness of spiky oscillations. The 3d volume 
is the region of spiky oscillations (Z > 2, defined in main 
text); colour gradient towards red shows increasing values 
of Z . A. The robustness for parameters A, B, S. The 2d 
plot shows the horizontal slice (white line) through the 3d 
volume at S = 0.029, the value used in Fig. The black dot 
corresponds to values of A, B used in Fig. [21 B. A similar 
plot for parameters A, B, C. 



The spikiness is extremely robust to variation of param- 
eters as shown in Fig. [5j Note the scales are logarithmic: 
parameters can be varied by several decades without going 
out of the oscillatory regime. The enormous robustness 
suggests that the spikiness of the oscillations is an impor- 
tant element in the design of NF-kB signalling and may 
be essential for its proper functioning as a transcription 
factor. We substantiate this idea in the subsequent sections. 

Sensitivity to IKK is high for spiky oscillations. 

Since IKK is the external signal to which the system 
responds, we begin by comparing the sensitivity of spiky 
and soft oscillations to changes in IKK concentration. We 
consider two quantities: the spike duration, defined as the 
amount of time N n spends above its average value, and the 
spike peak, defined as the maximum nuclear NF-kB con- 
centration during each cycle of oscillations. Fig. [HJA shows 
how the spike duration depends on IKK concentration. 
The sensitivity of the spike duration is very high in certain 
regions of spiky oscillations. It is especially large near the 
transition to soft oscillations. A similar sensitivity is seen 
in the peak NF-kB concentration (Fig. EJ3). Thus, the spike 
duration and peak are much more sensitive to (and therefore 
easier to regulate by) IKK for spiky than for soft oscillations. 
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Figure 6: Sensitivity to IKK. A. Spike duration, the frac- 
tion of time N n spends above its mean value, as a function 
of IKK concentration. The black dot shows the IKK value 
used in Fig. Blue and red signify, respectively, regions of 
spiky and soft oscillations. Notice the sharp response just 
before the transition to soft oscillations. B. Spike peak, the 
maximum concentration of nuclear NF-kB, as a function of 
IKK concentration. The black dot shows the IKK value 
used in Fig. |2 Blue and red signify, respectively, regions of 
spiky and soft oscillations. Notice the sharp response just 
before the transition to soft oscillations. 

Large Hill coefficients and regulation of down- 
stream genes. It is possible for genes regulated by NF-kB 
to inherit this sensitivity in the form of a high effective Hill 
coefficient. Consider a gene which has an operator site at 
which NF-kB dimers can bind and activate the gene: 

G + 2N ^ G*. 

koff 

To begin with, we assume that the binding of NF-kB to the 
operator is in equilibrium, i.e., k on and k a ff are much larger 




Figure 7: Equilibrium binding of NF-kB to a downstream 
gene. The plot shows the peak gene activity as a function of 
IKK concentration (open circles). The data has been fitted 
by a sigmoidal function of the form 1/(1 + (x/xo) h ). The 
least squares fit (solid line) gives an effective Hill coefficient 
h = 23.4. Inset: effective Hill coefficient obtained by similar 
fitting for different k ft/k on ratios. 

than the rates of all other processes in the NF-kB system. 
In that case the gene activity, G*, will follow the NF-kB 
concentration: 

N 2 

koff/kon +N%' 

Fig. shows the peak gene activity as a function of IKK 
concentration. The effective Hill coefficient of this response 
curve is over 20, much larger than the values obtained by 
typical ways of introducing cooperativity in gene regulation 
[241 125| . As the inset shows, the effective Hill coefficient 
remains above 20 for a large range of values of the ratio 
koff/kon i i- e -j genes controlled in this way show a very high 
sensitivity to the input signal. 

When k on and k a f f become comparable to other rates in 
the system the binding of NF-kB to the operator remains 
out of equilibrium. When fc // is small enough, the gene 
activity does not have enough time to decay completely be- 
tween spikes of NF-kB. Fig. [S] shows the peak activity as 
a function of IKK. In contrast to the equilibrium case, here 
the response is linear at best. Note also that the peak gene 
activity increases with IKK in contrast to the equilibirum 
case where it decreases. Thus, the same oscillations are ca- 
pable of regulating genes very differently, depending on their 
k n, of f values which are determined by their operator sites. 

In the extreme case where fc // is negligible, each spike 
of nuclear NF-kB contributes to increase the gene activity 
until it saturates to unity as shown by the solid line in Fig. 
If this gene controls the activity of a second one in a 
cascade: 

G* G* + p, 

G 2 + p'P G* 2 , 

then the latter will turn on later (dashed line in Fig. 0) with 
a time delay that depends on the timescales of transcription, 



Figure 8: Noncquilibrium binding of NF-kB to a down- 
stream gene. The plot shows the peak gene activity as a 
function of IKK concentration, with fc // = 0.001mm -1 
kept fixed while k on = 0.01/zM _1 min _1 for filled circles 
and 0.001fiAI~ 1 min~ 1 for open circles. 
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Figure 9: Strongly noncquilibrium binding of NF-kB to 
a downstream gene. The plot shows nuclear NF-kB con- 
centration (dotted line), the activity of a gene, G* (solid 
line), directly activated by NF-kB dimers with k on = 
0.1/uM~ 2 TOm -1 , k a ff = 0, and the activity of a second gene, 
G2 (dashed line), that is activated by the protein product 
of G (k t = O.Olmin , 7 P = lmin , see main text). 

translation and promoter activation. This is reminiscent 
of the experiments of rcf. |3j which show the gene IP-10 
turning on quickly after the introduction of IKK, while the 
gene RANTES turns on after a delay. 

Discussion 

What functional role, if any, do the oscillations in the NF-kB 
system play? There have been several suggestions: down- 
stream gene networks are perhaps regulated by the fre- 
quency of the oscillations, or the oscillations could be a 
by-product of rapid attenuation of NF-kB, or they might 
be used to make multiple evaluations of the input signal 
[261 177| . Barken et al. warn against overemphasising 
the physiological role of oscillations. Our approach to tack- 
ling this question has been to construct a reduced 3- variable 
model which, despite its simplicity, captures many charac- 



teristic features of the system. The simplicity of the model 
allows us to fully explore and understand the range of dy- 
namical behaviour it exhibits. In particular, we have shown 
that it is capable of both spiky and soft oscillations, and 
that the spiky oscillations are extremely robust to variation 
of parameters. 

The activity of genes downstream of NF-kB depends 
upon the amount of time for which NF-kB is present 
inside the nucleus in sufficiently large concentrations to 
dimerize and bind to those genes' operator sites. It seems 
reasonable to assume that NF-kB could signal different 
downstream genes, simply by regulating the amount and 
exposure time to IKK, provided the signalling system is 
sufficiently sensitive to changes in IKK concentration. We 
have shown that the spiky oscillations can, indeed, show 
a high sensitivity to IKK. This sensitivity allows a great 
versatility in the regulation of downstream genes by NF-kB. 
Where the cell requires a gene to be very sensitive to the 
IKK concentration, the NF-kB system can result in steep 
response curves with Hill coefficients larger than 20. And 
where a slower response is necessary, it can be achieved by 
adjusting the binding and dissociation constants of NF-kB 
to that operator site. Further, we found that cascades 
of different length could be used to turn certain genes 
on earlier or later. Given this versatility in regulatory 
strategics it seems likely that cells would have evolved to 
make use of these properties of the NF-kB oscillations. It 
remains for future experiments to uncover the particular 
ways NF-kB regulates specific genes. 

We thank J. Fcrkinghoff-Borg, E. Siggia and G. Tiana for 
useful discussions. This work was supported by the Danish 
Research Foundation. 
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Supporting Text 



7-variable model of NF-kB signalling 

We use the following abbreviations: N n !kN, free nuclear and cytoplasmic NF-kB; I m , IkB mRNA; free nuclear and 

cytoplasmic IkB; (NI) n &c(NI), nuclear and cytoplasmic NF-kB-IkB complex; IKK, IkB kinase. The 7-variable model is 
defined by the equations: 

—rjjr- = k m n N - k fn N n I n + k bn (NI) n , 

— = k t N n - 7 m i m , 

^ = hdm - k/NI + k b (NI) - kn n I + k Iout I n , 
dN 

— = -k f NI + (k b + a)(NI) - k Nin N, 
= k f NI - (k b + a)(NI) + k NIout (NI) n> 
= kjj n I — kj out I n — kf n N n I n + kb n {NI) n , 

^ ^ — = k fn N n I n - (k bn + k NIout )(NI) n . 

Figure IS1I shows a plot of nuclear NF-kB concentration obtained by integrating these equations using the follow- 
ing parameter values (Hoffmann ct al. (2002) Science 298, 1241): kmn — 5.4 min - , kun = 0.018 min~ 1 ,kj out = 
0.012 mm ,kj\riout ~ 0.83 min _1 ,/c t = 1.03 fiM~ l mm ,kti = 0.24 min _1 ,fc/ = kf n = 30 (J,M ~~ 1 min _ , kb = fc&n = 
0.03 min -1 , a = 1.05 x IKK min" 1 , j m = 0.017min -1 . The initial conditions were N = lfiM, IKK = 0.5/xM and all other 
concentrations zero. 



Reduction from 7-variables to 3-variables 

First, taking note of the fact that kf and kf n are large, we assume that all complexes are in equilibrium, i.e.: 

k f NI^(k b + a)(NI), 

k fn N n I n ps (kbn + k NIout )(NI) n . 

Simulations show that these are good approximations. In terms of 7*°* = 7„ + (NI) n and 7V* ot = N + (NI) = N tot — N n , 
which arc slowly varying, we can rewrite the above equations as follows: 

7 



(NI) = (N tot ~ N n ) 
N = (Ntot - N n ) 



Ki+r 

Ki 



Ki+r 



(NI) n =lf ot N ' 



In = L 



K N + N n 

tot K N 

71 K N + N n ' 



where Ki = (kb + ct)/kf = 0.035^M and = (fc&n + kpfj out )/kf n = 0.029/iA7, using the parameter values above. 
Using these expressions, the equations of the 7-variable model reduce to the following four (Fig. ISlf) : 



dN n (Ntot -N n ) I%*N, 

= ItNinKl — — KNIout 



dt " m 1 Kj + I iyiom K N + N n ' 
^In = k t N 2 - -y 7 

I, rX/ t iy n lra ± m-i 

dl (Ntot - N n )I 1% 

TT = k tl-l,n - Ct — — Kli„l + K Iout K N - 



dt Ki + 1 Km + N n 



7 variable model] 
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Figure SI: Schematic showing reduction of the 7- variable to a 3-variable model, as well as plots of oscillatory nuclear 
NF-kB from simulations of the 7 and 3-variable models using parameter values given in the text. 

JTtot Ttot Ttot AT 

ul n 1 T 1 n i 1 n iv « 

— Klin± — KioutlS-N— — — KNIout 



dt "K N + N n "'^K N + N n 

First, we note that the terms —kn n I and ki out Kjy K n +N in the dl/dt equation are much smaller than — a- — to ^ +J 

and can be neglected as long as IKK is non-zero. Secondly, simulations reveal that the term fcjv/o«t jy" + jv ; in the dl 1 ^ 1 /dt 
equation, also shows sharp spikes as a function of time which coincide with the spikes of N n . The value of this term is 
substantial only when N„ 3> K jy, i.e., during the spikes of 7V„, and at those times I 1 " 1 dips to its minimum. We therefore 
make the approximation that 7*°* can be replaced by its minimum value, /^ inl which satisfies the equation: 

Ttot Ttot AT 

1 T 1 jy n,min , n,min n 

Klini = Ki ou tK]y— + KNIouf 



In the regime where N n 3> K n this gives 



Using this we can reduce to a 3-variable model: 



K N + N n Kn + N„ 



rtot 
n.min 



*NIout 



dN n (N tot -N n ) IN n 

KNinKl — — Kir 



dt Ki + I S + N n 

dim , AT 2 T 

— = k t N n - 7 m 7 m , 
dl , T (Ntot - N n )I 

-77 = ktUrn - » • 

at Kr + I 



Rescaling the 3-variable model 

For ease of analysis, we reduce the number of parameters in the model by rescaling all variables to become dimcnsionlcss. 
We use the following transformations on the above equations: 



* - (l/7m)t, 



which gives: 



with 



N n ^N tot N n , 
I m -> (k t N? ot /j m )I m , 

i -> (k t k tl K t /^)i, 



A (l-N n ) B IN 



dt 



e + I 



S + NJ 



dt 



dl _ ^ (l-N n )I 
C e + I ' 



a = fc f7 A \ 7m « 0.007, 

kiinhhlNtot nr , _ 
±> = 5 w 954.5, 

/ 77 / 

«7m 



c 



hhiNtot 

Ntot 
Kill, ^ 



« 0.035, 
0.029, 



2 x 10 



-5 



Steady state solution of the 3-variable model 

The steady state values of N n , I m and I are solutions to 

A (l~N n ) B INn 
e + I 5 + N n 



K - I m = 0, 

i m -c^l = o. 

e + I 



I m and I can be eliminated using 



I 



Im = K, 



N 2 e 

± *n c 

C - CN n - N 2 • 



From this we find that the steady state value of N n is a solution of the equation 

BCe 2 



(C - CN n - N 2 ) 2 = 



A 5 + N n 



or equivalcntly, 

N* + (8 + 2C)N* + C 



2(5-l) + C-^ 



Nl + C[(C - 2)8 - 2C]N* + C 2 {1 - 2S)N n + C 2 8 = 0. 



In general, this has two real solutions, one with C — CN n — N 2 > and the other with C — CN n — N 2 < 0. The latter 
results in a negative value for / and therefore is not an acceptable solution. Thus we are left with only one fixed point. 



Linear stability of the fixed point 

We linearize the equations around the fixed point, which gives the Hessian 

n A(l-N n ) 





-1 

This matrix can be used to examine the stability of the fixed point . If Aj are the (possibly complex) eigenvalues of this 
matrix, then the fixed point is unstable if 

maxi [Re(Aj)] > 

and stable if 

max, [Re(Ai)] < 0. 

Linear production of I m 

In this section we consider the effect of taking the production of I m to be linear in N n , instead of N%: 

dN n {N tot -N n ) , IN n 

- k Nin Ki 



dt " m 1 Kj + I lin 5 + N n ' 
at 

dl , _ (N tot -N n )I 

Tt- ktlIm - a kT+i ■ 

Using very similar transformations, we rescale the variables: 

t (l/7 m )t, 
N n N tot N n , 

Im -> (ktN to t/lm)Im, 

I -» (k t k tl N tot /^ n )I, 



which gives: 



dN n = A (l-N n ) _ B IN n 



with 



dt e + I S + N n 7 

^=N -I 
at 

dl _ ^ (l~N n )I 
~df ~ m C e + I ' 

kNinKl"/ m 



A 



B = 



hhiNtot 

klinkt kfl 



/yd 

I n i 



c _«% 



5 = 



hhi ' 

K N 
Not 

c _ Kill 
kt ktiN tot 

If we use the same parameter values as before we do not get sustained oscillations. Not surprisingly, the region in 
parameter space where we get sustained oscillations has shifted. Simply taking A = 0.001, we get spiky oscillations, as in 
the original model (Fig. [S2l\ s |. Fig. IS2fe shows the response to changes in IKK. 



A. 



B. 




Figure S2: A. Sustained spiky oscillations of nuclear NF-kB when the production of I m is linearly dependent on N n . 
Parameter values: A = 0.001, B = 954.5, C = 0.035, S = 0.029 and e = 2 x 10~ 5 . B. Spike peak, the maximum 
concentration of nuclear NF-kB, as a function of IKK concentration. Blue and red signify, respectively, regions of spiky 
and soft oscillations. 




Figure S3: A. Steady state solution for / as a function of the parameter e, when the production of I m is linearly dependent 
on N n . Dashed line shows the region where the fixed point is unstable and the solid line shows the region where it is 
stable. Black dots mark the parameter values used in (B) and (C). The crossover from unstable to stable occurs where e 
becomes comparable to the steady state value of /, i.e., where the degradation of I stops being saturated. B. Trajectory 
of the system in the I-N n phase plane, for e = 1 x 10 -4 , which converges to a stable limit cycle. C. Trajectory of the 
system, for e = 2 x 10~ 4 , which converges to a stable fixed point. 



The steady state value of N n is now a solution to 

BCe 2 Nl 



(C - CN n - N n ) 



2 



A 8 + N n 

The steady value of I can be calculated from the value of N n using: 



I 



C - CN n - N n 



In general, there are again two real solutions, one with C — CN n — N n > and the other with C — CN n — N n < 0. The 
latter results in a negative value for I and therefore is not an acceptable solution. Thus we are left with only one fixed 
point. 

Linearizing the equations around the fixed point gives almost the same Hessian. The only difference is in one matrix 
element: J24 = 1 instead of 2N n . Fig. IS3I shows the stability of the fixed point as a function of e. 



